Collective T- and P- Odd Electromagnetic Moments 
in Nuclei with Octupole Deformations 



N. Auerbach/^'' V.V. Flambaum/^^ and V. Spevak*^^^ 
School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
^'^^ School of Physics, University of New South Wales, Sydney, 2052, NSW, Australia 
^ ; (January 28, 1996) 



O 



Abstract 
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Parity and time invariance violating forces produce collective P- and T- odd 
\ moments in nuclei with static octupole deformation. Collective Schiff mo- 

, ment, electric octupole and dipole and also magnetic quadrupole appear due 



to the mixing of rotational levels of opposite parity and can exceed single- 



. particle moments by more than a factor of 100. This enhancement is due 

Q"^ . to two factors, the collective nature of the intrinsic moments and the small 

energy separation between members of parity doublets. The above moments 
induce T- and P- odd effects in atoms and molecules. Experiments with such 
^ ' systems may improve substantially the limits on time reversal violation. 

> 

^ ' Parity and time invariance nonconserving nuclear moments induced by P-, T- odd nuclear 

■ forces were discussed e.g. in Refs. These moments can be enhanced in nuclei which 

have close to the ground state (g.s.) levels of the same spin as the g.s. but opposite parity 
0,^. An interesting possibility to enhance the effect is to consider mechanisms producing 
collective T-, P- odd moments. In Ref. [^] it was shown that the "spin hedgehog" mechanism 
produces a collective magnetic quadrupole. In the present paper we want to consider a 
different mechanism: mixing of opposite parity rotational levels (parity doublets) by T-, P- 
odd interaction in the nuclei with octupole deformation. This deformation was demonstrated 
to exist in nuclei from the Ra-Th and Ba-Sm region and produces such effects as parity 
doublets, large dipole and octupole moments in the intrinsic frame of reference and enhanced 
El and E3 transitions (see review [Q). 

Let us start our consideration from the expression for the electrostatic potential of a 
nucleus screened by the electrons of the atom. If we consider only the dipole T-, P- odd 
part of screening (Purcell-Ramsey-Schiff theorem [^]) one finds [Q: 
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«''('•) J3„ , 1 ,JT7^ r P-i'') ,,3. 



Here p(r) is the nuclear charge density / p{r)d?'r = Z, Ps{r) is spherically symmetric part 
of p{r), and d = J erp{r)d^r is the electric dipole moment (EDM) of the nucleus. The 
multipole expansion of (p{R) contains both T-, P- even and T-, P- odd terms. The dipole 
part in Eq. (|l]) is canceled out by the second term in this equation: 

- / e{rV^)p{r)d'r + \{dV)^ j p,{r)d'r = . (2) 

The next term is the quadrupole which is T-, P- even thus the first non zero T-, P- odd 
term is: 



If i 

(/ = --j ep(r)r„r^r^dVV„V/3V^ — 

+ ^('^V)V,V^-^ / Psi'^yo^rpd^'r . (3) 

Here raVpr^ is a reducible tensor. After separation of the trace there will be terms which 
will contain a vector S (Schiff) and a rank 3 Qap-y (octupole) moments |Q: 

3 3 3 

V = V'Schiff + V^octupolc 5 

where 

^ = ]^ (/ ^/^(^)^'^^'^ - / Ps{r)r'd'r^ (5) 
is the Schiff moment (SM) and 

QaM =J ep{r)[rarpry - ^{Sa^r^ + Sf^^r^ + 6a-yrp)]d^r 
2 2 [Att r 

Q... ^-Q, = -sIy j epiryYsod'r (6) 

is a tensor octupole moment. (Corrections to the Schiff and octupole field that arise from 
the fact that we use ps{r) instead of the full p{r) in Eqs. are of third order in the 

nuclear deformation). 

Here we will consider the collective SM, collective octupole and also collective dipole as 
well as the collective magnetic quadrupole resulting from the rotation of the dipole. The 
mechanism for collective SM, dipole and octupole is the following: collective moments in the 
body-fixed system of the deformed nucleus are assumed to exist without any T-, P- violation. 
However without T-, P- violation the average value of these moments for a rotational state 
in the laboratory system is zero. T-, P- odd mixing of rotational doublet states produces 
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an average orientation of the nuclear axis n along nuclear spin I. In the case of a nearly 
degenerate rotational doublet 

^^ = ^{\IMK)±\IM-K)) , (7) 

where K = In. If the T-, P- interaction mixes the members of the doublet with the 
coefficient a the total wave function is \t' = ^+ + or 

^ = ^{{l + a)\IMK) + {l-a)\IM-K)) , (8) 

one obtains 

, , , , KM 

{^\n^\^) = 2«^y:p^ . (9) 

The intrinsic electric dipole and Schiff moments are directed along n: d = dn and S = Sn 
and therefore have nonzero average values in the g.s. M = K = I. 

To elucidate the origin of collective T-, P- odd moments consider a simple "molecu- 
lar" model, shown in Fig. 1: two charges Zq and q with masses Z'm and m have coordi- 
nates Xi = —a and X2 = Z'a so the center of mass is at a; = 0. This "molecule" has dipole, 
quadrupole, octupole and Schiff moments. The electric dipole in this case is d = {Z — Z')aq. 
Consider Z' = Z, in this case d = 0, but the SM and octupole moment are not zero. The oc- 
tupole moment Qs is proportional to Z — 1, (for Z = 1 there is only quadrupole deformation). 
SM in the body-fixed frame is 

S = ^qa'Z{Z' - 1) . (10) 

For I = \ the SM 5* 7^ as opposed to the octupole moment which is 7^ in the 
intrinsic frame, but vanishes in the laboratory because one cannot satisfy angular momentum 
coupling. Thus the octupole deformation is hidden. It is possible to have a situation in which 
(i = 0, Q3 = in the laboratory system but the SM in the laboratory is S* 7^ 0. This result 
applies to any system, for example to an elementary particle (neutron, electron). Indeed, 
for spin s = | there is only one T-, P- odd formfactor [10,11|. However, we have shown that 



the two moments - EDM and SM are not necessarily related to each other. There is no 
contradiction here. The relativistic expression for the T-, P- odd electromagnetic current 
for s = I in momentum representation is 

= f{(f )i'ib(^^iviquij , (11) 

where q is the momentum transfer, 75 and a^^^ are Dirac matrices. The formfactor can be 
expanded 

f{q')=d + rm' + ... , (12) 
where d = /(O) is the electric dipole moment of the particle ^j. In the nonrelativistic limit 
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Jo = -f{q'W<Tq^ . (13) 
The electric potential produced by this current is 

^(g')=Jo^o. = -4vrz^(t/ + /W + -) • (14) 

where D^y = {A'Kg^y)/q^ is the photon propagator. In the coordinate representation 

,1-.^ ,1^5(r.), 

^(r) = d(TV^ + 47r/'(0)<TV5(r) + ... . (15) 

The first term in (y9(r) gives the long-range dipole field while the second term is the contact 
field of the SM, i.e. S ~ /'(0)cr (see Eq. (^)). Thus, the SM emerges from the same form- 
factor as the electric dipole. (Note the difference between the SM and the P-odd, T-even 
anapole moment which also produces a contact field but corresponds to an independent 
formfactor |TT]]). One can therefore have a priori a situation in which T-, P- is violated, the 



Schiff moment is not zero but the dipole moment of the particle is zero. 

The mechanism of rotational level mixing can also produce a magnetic quadrupole. In- 
deed, in the intrinsic frame of reference a deformed nucleus can have both a magnetic dipole 
and magnetic quadrupole without T-, P- violation. Then T-, P- odd interaction mixes ro- 
tational parity doublets and can produce magnetic quadrupole in the mixed state. It is also 
worth noting that higher T-, P- odd moments can appear due to rotation of lower moments. 
For example, rotating electric dipole produces magnetic quadrupole. However all these con- 
tributions to higher moments will be proportional to Lz/ (M^c) where is a large mass of 
the nucleus and consequently very small. 

The intrinsic moments of heavy deformed nuclei are well described using the two-fiuid 
liquid drop model |]12|-Q. We consider here even-odd nuclei, so electric moments, except 



the dipole, are determined by the moments of the even Z core. The surface of a deformed 
nucleus is 

i? = i?o(l + EA>«o) . (16) 

1=1 

The j3i deformation is determined from the requirement that the center of mass fixed at 
2; = 0, i.e. / zS'r = 0: 

V4Tfc5^(2( + l)(2i + 3) 
The proton density in case of deformed nucleus is fl 

Po Po e^Z 



P 



2 8 Ci?o 

' (18) 



X 



2 2 ^ Rq 2/ + 1 ^Ro' 
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where po = 3A/(47ri?Q) and C is the volume symmetry-energy coefficient of the hquid-drop 
model. The dipole moment generated by this proton distribution is in the lowest order of 
deformation [0-|T5] 



d 



intr 



eAZ- 



E 



(/2-l)(8/ + 9) 



•i+i 



(19) 



C 407r;e^ [(2/ + l)(2/ + 3)]3/2' 

The inclusion of the neutron skin effect as well as the shell correction reduces somewhat (iintr, 
nevertheless Eq. (^) with the C ~ 20-30 MeV fits experimental values quite well |12|,|14 



This moment appears only because the Coulomb force produces a relative shift of protons 
versus neutrons. The constant part of the density in Eq. ([ISD does not contribute to ciintr- 
The intrinsic Schiff moment turns out to be 



eARl 



'407r 



e^Z 19 
^70 



E 

1=2 



{I + i)AA+i 

(2/ + l)(2/ + 3) 



(20) 



Here the constant part of p in Eq. (|18|) gives the main contribution (about 90% in nuclei 
with Z ~ 90). The expression for the intrinsic octupole moment is [|n 



Q 



3 intr 



ezm 



■2V77r^' " 3 
The P- and T- odd potential has the form |§,16 



(/33 + ?t/^/?2/53 + ...) . 



(21) 



V 



PT 



^^poE-.(v./(..)) 



(22) 



where G = 10^^/m^ is the Fermi constant and Pt{r) = pof{r) is the nuclear density. 

We use here the particle-core model for a reflection-asymmetric nucleus [0,|1^. The T-, 
P- odd as well as P- odd T- even mixing was studied in this model recently ||19|. The wave 
functions in the model are [PTIM: 



ip 

MK 



2/+ 111/2 



167r2 



4(7r))DW<f^ 



(23) 



where R2{t^) denotes rotation through an angle vr about the intrinsic 2 axis. The $^ = $^ 
are particle-core intrinsic states of good parity p. Denoting the good parity core states x'^ 
and particle states 0'" we write 



(24) 



The states x'" are projections of the reflection-asymmetric states Xa EI 



(25) 
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The matrix elements of V^'^ are 



= a+6_(0+ + a_6+(0+ . (26) 

Note that the pseudoscalar operator V^^ cannot connect states of an even-even axially 
symmetric core Hill- The expectation value of a T-, P- odd operator O in a T-, P- 
admixed state is 

(l>+|6|<l+) = 2a,,(<|.+ |6|<|.r) + 25:a,,($+|d|$7) . (27) 

The matrix elements between core states are 

{X^\0\x-) = ixMx.,) ■ (28) 

Writing the one-body operator O as sum of core and particle parts O = Ocore + Op one 
obtains 

($+161$-) = (x^laorclXA)K i«- J + K ih- j) 

+ (0+|6p|0-)(a+ , + a_ ,6+ ,) (29) 

The contribution of the single neutron is small for the SM operator and is absent for 
the octupole moment. In the case of closely spaced doublets a+ j ~ a_ j, 6+ j ^ 6_ j and 
(a+ ja_ j + 6+ j6„ j) ~ 5ij. The expressions for the expectation values of a T-, P- odd op- 
erator of rank / in the body-fixed and laboratory (for the g.s. I = M = K []T3[) systems 
becomes 

{^itK\o\^itK) = {im\n)'{^t\o\^t) ■ m 

Currently, the best limits on SchifF moments and the coupling constants of T-, P- violat- 
ing nucleon-nucleon interactions are obtained from the measurements of the electric dipole 
moments in ^^^Hg, ^^^Xe I^Uj atoms and TIF molecule [PT| , P^ . Nuclei of these atoms do not 
have octupole deformation. However, similar experiments can be done with heavy atoms 
(Ra, Rn) which are electronic structure analogs of these atoms but their nuclei have octupole 
deformation. 

Our calculations were performed for relatively long lived even-odd isotopes ^^^'^^^Ra and 
^^^Rn. Variants of the model used here are shown to describe quite well the g.s. parity 



doublets in the Ra-Th region ||T7| , |23| j8[] . We used here the same version as in Ref. [|l9 
the deformation and core parity splitting parameters were taken from Ref. [|T^,^. The 
calculations of the mixing coefficients were performed using Nilsson potential. The ^^^'^^^Ra 
and ^^^Rn have the g.s. |^ and |. The octupole deformation for all these isotopes is 
~ 0.1. Our results are shown in table |I[ The mixing coefficients in our calculations 
are in the range (0.6-7.) xlO^''?]. The Schiff moments of the reflection asymmetric nuclei 
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in the intrinsic system are in the range (22.-29.)efm^ Correspondingly the T-, P- odd 
Schiff moments (Eq. (pO[)) are in the range (2.-20.) xlO~^r7efm'^. This is about two orders of 
magnitude larger than the largest single-particle estimate of 4.xlO~^r7efm^ given in Ref. 0] 

for 237Np. 

Since the electronic structure of Ra is similar to that of Hg and Rn is similar to Xe, we 
can use results of atomic calculations for Hg, Xe |p^,|5[| and write 



4t(Ra) = 4t(Hg) 



{SZ^R 



1/2 jRa 



{SZmi/2h, 



4.(Rn) = 4t(Xe)i§;m^ , (31) 

where we have taken into account the Z- dependence of the Schiff moment contribution 
to the atomic EDM (iat due to increase of electronic wave functions near the nucleus |^ 
(iat ~ SZ^Ri/2. The relativistic factor is given by 

U,,. ,^J''^''J '-EIky-''-\ (32) 



|r(2Tv2 + 1)1= V 



a 



B 



where 71/2 = [1 — {ZaYY^"^ , the Bohr radius and Rq is the nuclear radius. 

We made an estimate also for ^^^Pa which has the smallest energy splitting between 
members of g.s. parity doublet among the known octupole deformed nuclei Q. We assumed 
the same atomic physics parameters as for Ra, Rn. As seen in table | the EDM of Ra 
and Rn are 10^-10^ times larger than EDM of Hg and Xe. We should stress that the phe- 
nomenon described in our work occurs probably in many other nuclei that possess octupole 
deformation (for example in the Ba-Sm region) and which might be more suitable choice 
for experimental studies. Experiments with atoms or molecules containing these nuclei may 
improve substantially the limits on time reversal violation. 
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TABLES 

TABLE L Admixture coefficients a (absolute values), experimental energy splitting between 
the g.s. doublet levels AE = E~ — , intrinsic Schiff moments and Schiff moments as well as 
induced atomic dipole moments. The values for ^^^Hg and ^^^Xe from Ref. [^J5| are given for 
comparison. 
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FIGURE 1. A "molecular" model of octupole deformation: two charges Zq and q with masses 
Z'm and m placed at —a and Z'a with respect to the center of mass. 
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